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Applying perturbative QCD methods recently seen to give a good description 

> ■ 

of the two body hadronic decays of the B meson, we address the question of 

y ■ 

bound-state effects on the decay B — > i^*7. Consistent with most analyses, we 
demonstrate that gluonic penguins, with photonic bremsstrahlung off a quark, 
change the decay rate by only a few percent. However, explicit off-shell b quark 
effects normally discarded are found to be large in amplitude, although in the 
standard model accidents of phase minimize the effect on the rate. Using an 
asymptotic distribution amplitude for the K* and just the standard model, we 
can obtain a branching ratio of a few xlO -5 , consistent with the observed rate. 
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1 Introduction 



This note reports on potential bound state modifications to the analysis of B — > K*^, 
which is usually given solely in terms the on-shell subprocess b — > 37. 

The flavor changing neutral currents involved in the decays of the B into K*^ do 
not exist to leading order in the standard model, but can occur in second order in the 
Weak interaction via emission and reabsorbtion of W bosons [I]. These loop diagrams 
are often called "penguins," and their magnitude can be greatly modified by strong 
interaction effects 0, 0, |J. 

There is recent further interest in these decays because additional penguin-like con- 
tributions could come from particles not in the minimal standard model 0. Contri- 
butions to B into K*^ decay from loops of non-standard-model particles (such as loops 
of supersymmetric particles that might be called "penguinos" ) would be a signal of their 
existence. To take advantage of this possibility, more precise study of the decay in the 
standard model needs to be undertaken. 

The subprocess b — > 57, taken as a free decay, is usually treated as the only fla- 
vor changing contribution leading to B — > K*j 0. However, bound state effects could 
seriously modify results coming from this assumption. Bound state effects include mod- 
ifications due to the quarks being off shell in b — ■> 57, contributions involving gluonic 
penguins or double (photon plus gluon) penguins, and contributions from annihilation 
diagrams. The latter involve no neutral flavor changing currents at all. 

We use shall perturbative QCD (pQCD) in our analysis (see also [§]), a methodology 
we have previously applied [0, [T0[] to hadronic decays and semileptonic form factors of 



the B, inspired by Ref. Examples of the Feynman diagrams we calculate are given 
in Fig. [I]. We require as input the effective vertices that result from the penguin diagram 
analyses 0, |J. Thereafter our calculations are quite explicit and are detailed below. 

A preview of our results is as follows: Diagrams involving the subprocess b — > 57 do 
dominate, and keeping just contributions from the most commonly cited effective vertex 
gives close to the correct answer. However, some luck underlies the last statement. Since 
the internal heavy quark propagator can go off shell, there is an additional, independent, 
effective vertex |3| that can contribute. It does so with an amplitude whose magnitude is 
about 2/3 that from the usual vertex, and only because of fortunate phase relations is the 
magnitude of the sum nearly the same as for the usual vertex alone. Other diagrams are 
shown to be small, although they lead to a someday measurable few percent difference 
in the decay rates of the charged and neutral B into K*j decay. 
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A diagram we are forced to omit for now is the double penguin, Fig |l]e, where both a 
photon and a gluon come out of the loop, as the effective vertex it gives is not calculated. 



2 Calculations 

We now begin to describe our calculations in more detail. The penguins are represented 
by blobs in Fig. [I], which may be interpreted as an effective Hamiltonian, expanded as 

Heff = -£fiV*VZ £ C^)O t (fi). (1) 

The operators Oj are listed in references || and ||. 

Consider the photonic penguin diagrams in Fig. |l|a. If the incoming and outgoing 
quarks in b — > sj are on shell, there is only one relevant operator in H e ff. Refs. Q 
contain two operators that can contribute, but using i]t)q = m q q they can be seen to be 
equivalent. We will write the operators so that only one contributes when the quarks are 
on-shell. The commonly used operator relevant to radiative b decay is then 

° 7 = 1^2 m ^ F ^\ (! + 7s) b. (2) 

The numbering is that of Ref. ||; unfortunately the notations of || and Q do not match. 
The other operator is 



where the numbering is from Ref. 0] and the prime reminds us that we have put the 
on shell part into the other operator (and that we changed the location of a factor 

Q d = -1/3). 

Let us here record the coefficients at the W mass scale, 
C' 2 {m w ) = — - ((1 - x)(18x 2 - liar - 1) + 2(3x - 2)(5x - 2) lnx) = -0.47, (4) 



24(1 - x) 4 

C 7 (m w ) = —-^ — - f(l - x)(8x 2 + 5x - 7) + 6x(3x - 2) Inx) = -0.19, (5) 
24(1 — x) 4 v ' 

where the numerical values are for x = m 2 /w^ = 4, and also record how these operators 
evolve down to lower scale, 

CM = C' 2 {mw) - | (l - V' 2/Po ) ~ f x {V 4/P ° ~ l) , (6) 
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CM = V- Wm [Cr(mw) - ^ (V° /3/3 ° - l) - ^ (v W ~ 1 



(7) 



Quantity r\ is a s (n) /a s {mw) and (3$ is 11 — (2/3)ny, and following standard practice we 
have neglected mixing with operators that give small effect. Both CV and C 2 increase in 
magnitude with decreasing scale \i. 

We make the peaking approximation for the distribution amplitude of the B 
meson, wherein 

<I>b{x) = -^=S{xx-€b). (8) 

The decay constant of the B is fs and x\ is the light cone momentum fraction carried 
by the light quark. The mass of the B is given by = + A# and e# = Ab/tub- For 
the K* distribution amplitude, we write 

<f>K*{y) = "/3fK*yiy2(f>K*(y)- (9) 

The normalization is 



dy 1 <t> K .{y) = ^ m (10) 
o 2V3 

so that 4>K*(y) is unity for the (super-)asymptotic distribution amplitude. We also make 
the approximation that m^* = 0. The second diagram of Fig. |l]a is then zero. 

The spin projection operators for the initial and final hadronic states are 75^ — 
itib)/\/2 and + (mj<-«))/\/2, respectively, with p, k, and q being the momenta of 
the B, K*, and photon, and £ the polarization vector of the K*. Angular momentum 
conservation allows only transverse polarizations. 

The result is 

M photonic penguin = — (C 7 (ji)I - 2CMI') (p-qe-£ + ie^ a ^q»e a e) , (11) 

rriB^B v J 

where e is the polarization of the photon and 

G = C F ^^=V tb V t :f B f K * (12) 

with C F = 4/3. Also 

I = (l-e B ) dy l( f) K *(y) — 13 

Jo Hi — 2e B — t0 + 



I ^)(-i + (l-2c B ) 



«7r + In 

2e B 



(14) 
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and 

I'= [ 1 dy 1 (l-y 1 )$K-(v) = L (15) 
Jo 2 

where the integrated results are for the asymptotic distribution amplitude. The imaginary 
part come from an internal propagator going on shell. This is often taken as a signal 
that pQCD is inapplicable. However, if properties of the reaction overall dictate short 



distance propagation only, then pQCD can still be used |12| . This is the situation here, 
as discussed in |J and 

Numerical results will be given after we have discussed what prove to be the sub- 
dominant contributions. Our understanding of why they are subdominant is helped by 
some order of magnitude estimates. The crucial expansion parameter is 1/sb and its log- 
arithms. Factors of e# come from the propagators, and can also be induced, depending 
on circumstances, by the factor yiy 2 in the K* distribution amplitude. 

In the photonic penguin diagrams, Fig. |I|a, the gluon connects at the lower vertex 
to on-shell quarks and its propagator gives a factor proportional to 1/yies, where y± 
and €b are the momentum fractions of the two light quarks. Thus appears one factor of 
l/es'i the yi is canceled from the K* distribution amplitude. The b quark propagator 
is involved in two subprocesses: scattering from the light quark by gluon exchange and 
decay into the s quark plus photon. Both are possible for an on-shell b quark, and the 
b quark does go on-shell in this diagram when yi = 2e#. The b quark propagator thus 
contributes an imaginary pole term and a real principal value term (or just a real term, 
for the operator 0' 2 ) to the integral involving the K* distribution amplitude, and one of 
them gives (roughly speaking) an in and the other gives a logarithm of l/eg. Now we 
have accounted for the e# factors that appear in the photonic penguins, 

M photonic penguin « (factors) x — x (c 7 (ji) x O (in or In— ^ + C' 2 {ji) x 0(1)) . (16) 

The gluonic penguin graph with emission of the photon from the b or s quark (Fig. 
|T|b) does not allow the quark propagator to be on shell. For example, in the diagram 
with photon emission from the s quark, the internal s quark decaying into an on-shell 
photon and and on-shell s quark must have a momentum squared larger than the mass 
squared of the s. A similar argument shows the b propagator is never on shell in the 
diagram with photon emission from the b quark. The gluon propagator is still spacelike 
and still gives a l/e#, but compared to the previous case we lose the log(l/e.e) or in that 
came from quark propagator, 

M h or s bremss w (similar factors) x — x C 8 (/i) x O (1) . (17) 
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As we shall see, another significant reduction comes from the replacement of coefficient 
Cf by its gluonic counterpart Cg. 

For the spectator bremsstrahlung case, Fig. p]c, the quark propagator cannot go on- 
shell. However, the gluon propagator can. A factor (l/e B ) that came from the gluon 
propagator is lost, and replaced by factors m or log(l/ee) that come from integrating 
the K* quarks's momentum fraction over the gluon pole, yielding 

M spectator bremss « (similar factors) x (e B )° x C 8 (/i) x O (in or • ( 18 ) 

Thus the gluonic penguin diagrams are suppressed by powers of or logs thereof, as 
well as by the ratio C 8 /Cj. 

For the actual calculations involving the gluonic penguin we kept just 8 in the 
effective hamiltonian, where 

°% = S—2 ^so^G^t}- (1 + 75 ) b. (19) 

The numbering is that of Ref. 0. Other operators are possible when the gluon or 
quarks are off-shell. We have not explicitly calculated their contributions in this case (in 
part because their evolution has not been calculated), but have verified that the order of 
magnitude estimates are not upset, i.e., they are not leading in l/e B . 
The results are 

M b or s bremss = ^^C B (ji)I' (p-qe-£ + i^ ua ^q v e a e) , (20) 
m B e B v ' 

where = —1/3 and G and I' have the same meanings as before, and 

M spectat or bremss = C % (ll)I Q (e B ) (p • qt • f + ie^a^ <f 6° ^) , (21) 

m B v ' 

where e q is the quark charge for the spectator. We neglected some small terms in the 
numerator and let 

I Q (e B ) = -m + hx^—^-. (22) 

The coefficient C 8 is 

C 8 (m w ) = ((x - l)(x 2 - 5x - 2) + 6x\nx) = -0.094, (23) 

8(1 — xp v ' 

where the numerical value is again for x = 4, and it scales like, 

(24) 



c7 8 (a0 = \ C8{mw) _ n ^ m _ ^ + 35 ^ m _ ^ 
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As the renormalization scale decreases, the magnitude of C$ decreases, actually passing 
through zero at \i ~ 6 GeV. 

Additionally, there are the annihilation graphs, Fig. |l|d. These can contribute only 
to decay. To leading order in (e#) -1 the result comes from bremsstrahlung off the 
initial u quark and gives 

(p-qe-Z + ie^p^e^), (25) 

where we kept rriK* when it appeared as an overall factor and where e u = 2/3. The 
quantity in square brackets differs from the quantity G used earlier in lacking the strong 
interaction factors CVa s /47r and in having different CKM factors. While it is interesting 
that the decay proceeds at all without flavor changing neutral currents, the result turns 
out small. Not having the gluon exchange is a plus numerically, but the factors mK*/fn B 
and V u b ensure the small result. 

3 Numerical results 

The numerical results should not depend on the renormalization scale. However, as we 
are most familiar with the wave functions or distribution amplitudes at a typical hadronic 
scale, say /i ~ 1 GeV, we should evaluate the other quantities at the same scale. There 
are very big changes in the Ci from their values at the W mass scale. 
For the sake of definiteness we shall use 

A QCD = 100 MeV, 

A B = 500 MeV, 
my/ = 81 GeV, 



m t = 


2m w , 


V ts = 


-0.045, 


v tb = 


0.999, 


v uh = 


0.0045, 


V us ~- 


= 0.22, 



tb = 1.46 picoseconds, 



M 9 



2e M m K * 
m B e B m B 



G F 
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f B = 132 MeV, 

and 

f K * = 151 MeV 

Our convention has f„ = 93 MeV and the signs of the CKM parameters follow a "stan- 
dard" advocated in [[14]]. The sign of V t bV t */V u bV* s is what is crucial and it does not 
depend on conventions. We also use the asymptotic form for the K* distribution ampli- 
tude and will mention results with another form later. We will express each contribution 
to the amplitude as 

Mi = U x ~ (e • £ + z(p ■ qyh^^e^e) , (26) 

whereupon 

1 W (27) 



167T7«b 

for t being the sum of the (and neglecting the K* mass). 

The scale we should use should be compatible with the scale that our wave functions 
and distribution amplitudes are determined at, and this is turn should be consistent with 
the scale of the four-momentum squared of the off-shell gluon. This suggests \i ~ 0.8 
GeV, which we shall use. We extrapolate the coefficients according to the formulas given 
earlier. Much of the change due to the extrapolation occurs as the scale changes from 
my/ down to order m_e, at least for the large terms Cj and C 2 . For the running coupling 
we use a s = 4tt / p \ia(fi 2 / Aq CD ) with number of flavors appropriate to the scale. 

From operators 7 and 0' 2 we get 

t 7 = _o.95 -3.56i 

t v = 2.40 (28) 

with the contribution from photonic penguins being the sum of these two. The amplitudes 
are in units of 10 -8 GeV. For others we get 



or s bremss 0.08, 
^spectator bremss 0.04 0.05i, 

tann = 0.06. (29) 

The last two are given for the charged B. The photonic penguins dominate, although 
the other graphs contribute circa 10% corrections to the real parts of the decay amplitude. 
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The extra operator that we have considered in the photonic penguin calculation gives 
an amplitude with a magnitude that is about 2/3 as large as the amplitude from the 
operator normally considered. It is the luck of the phases that keeps its effect small. It 
increases the rate by not even 10% since roughly speaking what it does is just change 
the sign of the real part of the amplitude. 
Specifically, we get a branching fraction 

Br(B -> K*j) = 1.24 x 1(T 5 (30) 

including just the photonic penguin terms. Keeping only the most usual operator Oj 
would reduce the branching fraction to 1.13 x 10 -5 . It seems inconsistent to include 
the smaller contributions since they may be smaller than the errors induced by our 
approximations upon the big terms. However, keeping all terms anyway gives 1.24 x 10~ 5 
(unchanged from above) for the neutral B and 1.31 x 10~ 5 for the charged B. The relative 
size of the neutral and charged B decays should be about right and would be interesting 
to observe as more precise data becomes available. 

It is possible that the distribution amplitude for the transversely polarized K* is 
narrower than the asymptotic one. Chernyak, Zhitnitsky, and Zhitnitsky [T5| suggested 
a distribution amplitude 

My) = Hym) 2 , (3i) 

albeit this was for the transverse p. If we use this distribution amplitude for the K*, our 
calculated branching fractions are roughly halved. 

The choice of Aqcd was made consistent with some of our own earlier work H 10], 



but could be varied (the earlier situation is much less sensitive to the value of this 
quantity than the present case will prove to be). If we let Aqcd = 200 MeV, leaving 
other parameter choices untouched, the branching ratio with the asymptotic distribution 
amplitude changes to 

Br(B -> K*j) = 3.50 x 10~ 5 (32) 

including just photonic penguins terms, with commensurate changes in results keeping 
all terms or just Of. These values are in accord with present experimental data |[L6|| . It 
is also clear that values for other parameters could be varied somewhat from values that 
we have used. 
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4 Conclusion 



It seems with present knowledge, the actual B — > i^*7 decay rate is sensitive to pa- 
rameters of the bound state and to parameters governing the evolution of QCD. Still, a 
number of conclusions may be drawn. 

Contributions from gluonic penguin and annihilation diagrams — which contribute to 
the physical decay but not to b — > 57 — have been calculated here. They change the decay 
rate by a few percent and so are not worrisome until the experiments are considerably 
more precise. 

Also calculated here, and more significant, are effects due to decay quarks being off- 
shell. This brings into play another operator in the effective hamiltonian for b — > sj, and 
this new operator produces an amplitude of noticeable magnitude. However, its phase is 
such that the effect on the decay is under 10%. 

Regarding the future, there is a need to calculate the double penguin diagrams men- 
tioned in the introduction and illustrated in Fig. |l|e, including the QCD corrections to 
them. For now, we can at least note that these diagrams will not contribute an imaginary 
part to the amplitude, so that the imaginary part that is already there puts a lower bound 
on the total result. Also, one will wish to eventually study the totality of B — > since 
this is closer to the basic b — > 37 than any individual exclusive channel. Still, the physical 
B — > Xg'-f always involves spectators and always has a contribution from operators that 
only contribute if the b is off-shell. We have not calculated for situations other than the 
lowest "A s " and do not know what size corrections ensue overall, or even if the phase 
situation for the off-shell contributions persists. 

None-the-less, the opportunity to test the flavor changing neutral currents induced 
in the standard model and to search for evidence of particles or phenomena beyond the 
standard model makes B decay into K* / y and into A s 7 interesting, and makes calculations 
to determine precisely the standard model contributions to these decays a worthwhile and 
necessary pursuit. 
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Figure caption 

Fig. [I]. Diagrams. The effective vertex, due for example to a W and t quark loop, is 
represented as an oval blob; (a) shows the photonic penguins, (b) shows gluonic pen- 
guins with bremsstrahlung from the b or s quark, (c) shows gluonic penguins with 
bremsstrahlung from the spectator quark, (d) shows two of four lowest order annihi- 
lation diagrams that could give charged B —>■ K*j, and (e) shows the "double penguin," 
once as a blob and once showing one example of a contribution to the blob. 
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Figure 1: 
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